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1. Introduction 

Stringent tests of the Standard Model and the search for new physics pass through a detailed 
study of physical processes involving the ft-quark. It is clear that, despite the present (B-factories, 
LHCb) and future (SuperB factories) experimental programs, lattice QCD results of hadronic pa- 
rameters need to have reduced uncertainties at the level of ~ 1 % (see ref s [1,2] for recent reviews). 
In this proceeding we report on the ongoing project, within the ETMC collaboration, to compute 
B-physics hadronic parameters with Wilson twisted mass fermions. We extend and improve our 
previous analysis [3, 4] in two ways: 

• we optimize the interpolating operators for heavy-light systems to better project onto the 
fundamental state 

• we extend the range of heavy masses, /f/,, considered reaching values of approximately /X/j ~ 
2.5 \l c , where \i c is the charm quark mass. 

2. Improved interpolating operators 

We use, for this analysis, the Nf = 2 dynamical gauge configurations with up and down mass degen- 
erate quarks, \i u j d = /I/, generated by the European Twisted Mass Collaboration (ETMC) [5]. The 
lattice action is the tree-level improved Symanzik gauge action [6] and the twisted mass quark ac- 
tion [7] at maximal twist [8]. The strange, jl s , and the charm (heavy), jx c (jU/,), quarks are quenched 
in this work. 

We have used four lattice spacings a = {0.098(3), 0.085(2), 0.067(2), 0.054(1)} fm [9] correspond- 
ing to jS = {3.80,3.90,4.05,4.20}. We also use the values of the renormalized quark masses 
V-u/d = 3.6(2) MeV and Jl s = 95(6) MeV [9] and the pseudoscalar density renormalization con- 
stants Z|F(2GeV) = {0.411(12), 0.437(7), 0.477(6), 0.501(20)} [10] at the four beta values. We 
denote by a "bar" the quark masses renormalized in the MS scheme at a renormalization scale of 2 
GeV. The values of the bare valence quark masses used in this calculation are given in tab. 1 . 

To keep the noise-to-signal ratio under control we extract meson masses at relatively small 
temporal separations. To improve the projection onto the fundamental state, we need trial states 
with large overlap with the lowest lying energy eigenstate. We have constructed a smeared inter- 
polating operator using Gaussian smearing [1 1], i.e. we have inverted on a source <I> S 

* s «(l + Kfca 2 V£ ra ) Aro <fr L , kg = 4 N G = 30, (2.1) 

where <I> L is the standard local source and Vape is the lattice covariant derivative with APE smeared 
gauge links (0!ape = 0-5, Nape = 20) [12]. We have used both smeared sources and sinks finding 
that the best correlator in respect of signal to noise ratio is the smeared-local (SL) correlator where 
the source is smeared and the sink is local. 

To further improve the overlap with the ground state of the B-system we have constructed the 
source 

<£(w)°c ft><J> s + (l-ff>)<J> L , (2.2) 

dependent on a tunable parameter CO. We have computed correlators where the source is <J>(o>) and 
the sink is smeared (ooS) for several values of ft). We then optimize the value of ft) to achieve a 
projection onto the ground state at earlier Euclidean times than with SL conelators 1 . In the left plot 

'Details on the optimization procedure will be given in a forthcoming publication. 
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a\l s 


a\x h 


3.80 


0.0080, 0.0110 


0.0175, 0.0194 
0.0213 


0.1982, 0.2331, 0.2742, 0.3225, 0.3793 0.4461 
0.5246 


3.90 


0.0030, 0.0040 
0.0064, 0.0085, 0.0100 


0.0159, 0.0177 
0.0195 


0.1828, 0.2150, 0.2529, 0.2974, 0.3498 
0.4114, 0.4839 


4.05 


0.0030, 0.0060 
0.0080 


0.0139, 0.0154 
0.0169 


0.1572, 0.1849, 0.2175, 0.2558, 0.3008 
0.3538, 0.4162 


4.20 


0.0020, 0.0065 


0.0116, 0.0129 
0.0142 


0.13315, 0.1566, ,0.1842, 0.2166, 0.2548 
0.2997, 0.3525 



Table 1: Bare parameters used in this computation. The volumes L? x T simulated are L/a 
<; 24. 24. 32. 48;. with T/L = 2, respectively for j8 = {3.80,3.90,4.05,4.20}. Only for the ensemble with 
a\li = 0.003 and j8 = 3.9 we have L/a = 32. 



of fig. 1 we show the Euclidean time dependence of the effective masses for a particular simulation 
point obtained from local-local (LL), SL and coS correlation functions. The plot describes well 
the general properties for heavy-light effective masses: the SL correlators are already substantially 
improved for a safe extraction of the ground state mass compared with the LL correlators. The coS 
correlation functions perform even better. This can be fully appreciated in the right plot of fig. 1 
where we compare for different plateau regions the effective masses extracted from the SL and the 
coS correlation functions. It is clear that with coS correlation functions we can extract the mass of 
the ground state using results at Euclidean times smaller than with the SL ones. In the same plot 
we also compare these two determinations with the one obtained with the GEVP [13], obtaining 
perfectly consistent results within statistical uncertainties. We prefer not to use the GEVP in our 
final analysis because we have found that it gives noisier estimates for the pseudoscalar decay 
constant. All the results presented in the following sections take advantage of these improvements. 

3. The £-quark mass and decay constants /g and fs s 

To determine the Z?-quark mass, fs, fg s and the bag parameters, we implement the ratio 
method [4]. We refer to this ref. for the details of the method. Here we briefly recall the basic 
steps for the determination of the Z>-quark mass. For the other quantities discussed in these pro- 
ceedings the same strategy applies. The method is suggested by the HQET asymptotic behavior of 
the heavy-light meson mass Mm in the pole heavy-quark mass fif° le 

We consider a sequence of heavy-quark masses Jl^}' < Jlf 4 ' < • • • < fv^' with fixed ratio A, i.e. 
~Ph^ = • The k ev quantities of the method are ratios of heavy-light meson masses at sub- 

sequent values of the heavy-quark mass, y(^ h properly normalized [4]. From eq. (3.1) 

and QCD asymptotic freedom it follows that the ratios y(~p h , A; ju"/,a) have an exact static limit: 

Jim limy(/l)"\A;/I,,a) = 1 . (3.2) 
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Figure 1: Left plot: Euclidean time dependence of the effective masses at /3 — 3.8, aji/ — 0.08, ajj,/, — 0.5246 
obtained from LL, SL and ©S correlation functions. Right plot: effective masses at the same simulation point 
as for the left plot extracted from the SL and the 0)S correlation for different plateau ranges. The ranges, in 
units of the lattice spacing a, are indicated in square brackets. We also plot for comparison the result of a 
GEVP analysis. 




The value of X is chosen in such a way that after a finite number of steps the heavy-light meson 
mass assumes the experimental value Mg = 5.279 GeV. In order to implement this condition, the 
lattice data at the four lattice spacings are interpolated at specific values of the heavy-quark mass. 
In the left plot of fig. 2 we show the chiral and continuum extrapolation of the ratio of the heavy- 
light meson masses evaluated at the heaviest quark mass. A phenomenological fit linear in the 
light-quark mass and in a 2 described the lattice data rather well. We observe that discretization 
eiTors are well under control. The values of the ratios y(p h ,X) extrapolated to the chiral and 
continuum limit have a non-perturbative heavy-quark mass dependence that is well described by 
the HQET-inspired function 

; fel) = l + f + (33) 

as it can be seen from the right plot in fig. 2. From this plot it is clear that the improvements we 
have implemented in this analysis allow to have an accurate description of the non-perturbative 
dependence on ~p h of the ratio y. In particular there is little doubt that the numerical data are 
well described by the fit ansatz in eq. (3.3). The formula (3.3), with the fit parameters T]i and T]2, 
provides the non-perturbative description of the heavy-quark mass dependence of the ratio y(]Z h ,X) 
over the whole range of heavy-masses. We can thus compute the value of the Z;-quark mass simply 
using the knowledge of y(p h ,X) to iterate the value of the triggering heavy-light meson mass 
Mhi(Ji h ) from the charm region to the experimental value of the B meson mass M#. This is 
always possible with a slight tuning of X and Wh • One finds X = 1.1784, Jv^' = 1.14 GeV => 
J[ b = X K Jl^ (K = 9), that leads after a 4-loop evolution to 

nF*(Hb)\ Nf =2 = 4.35(12) GeV. (3.4) 

This result is consistent with our previous determination [3] with slightly reduced error, and a fully 
equivalent one is obtained if we use as input the heavy-strange meson mass Mf ls . The 2.7% relative 
eiTor is dominated by the ones associated to the scale and Zp determinations. Note that a Z?-quark 
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Figure 2: Left plot: chiral-continuum extrapolation of the ratio of heavy-light meson masses at the heav- 

in Q\ 

iest quark masses, jj, h ' . Right plot: non-perturbative heavy-quark mass dependence of the ratio function 
y(]J h ,X), for X = 1.1784 in the continuum. For comparison we show our previous results [3] without the 
interpolating operators improvement discussed in sect. 2. 



mass value smaller by ~ 120 MeV would be found by setting Nf = 4 (rather than 2) in the evolution 
from 2 GeV to the Z?-scale. The statistical and systematic uncertainties coming from the application 
of the ratio method turn out to be negligible, after the improvements just discussed. 

To compute heavy-light and heavy-strange decay constants a completely analogous strategy 
can be adopted. We find advantageous to define ratios with exactly known static limit for the heavy- 
strange decay constant, z s (jl^, A.;jU;,jU 4 .,a), and for the ratio /bJ/b, C(Ph i^'->V-hT^si a ) P]. As 
for the fr-quark mass, fhs at the triggering point and the ratio z s have a smooth chiral-continuum 
extrapolation with cutoff effects always well under control. The heavy-quark mass dependence in 
the continuum of z s can again be described by a formula as the one in eq. (3.3), as shown in the left 
plot of fig. 3. The application of the ratio method using as input the value we determined for JI b 
leads us to 

/fij = 234(6) MeV. (3.5) 

The chiral-continuum extrapolation for the double ratio £(jLj" , X;Jii,Jl s ,a) does not pose any 
problem and its heavy-quark mass dependence is very weak easing the interpolation at the ft-quark 
mass. At the triggering mass, p.j^\ the data for fhs/ f hi as a function of p.i suggest that the /j-quark 
still behaves as a relativistic one. Thus for the chiral-continuum extrapolation of fhs/ fu(P-jp) we 
find convenient to study the double ratio [fhs/fhi] • [fsi/f^ P ] , for which a smooth dependence on 
the light-quark mass is expected [14, 15]. We find that a linear fit in the light-quark mass describes 
well the data as it can be seen from the right plot in fig. 3. This leads to the values 

^ = 1.19(5), fs = 197(10) MeV, (3.6) 

JB 

where we add in quadrature the statistical and systemtic uncertainties. 

4. Bag parameters 

The B-parameter of the renormalized operator in QCD is defined as 

(5,|Of= 2 |S ? r^|/isl ? (MK ^ 
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Figure 3: Left plot: same as the right plot in fig. 2 fot the ratio function z s {li- h ,X). Right plot: chiral 
extrapolation for the double ratio [fhs/fhi] ■ [fsi/f^] at the triggering point, ~p h . 



where 0^ =1 = [&7^(1 — Ys)q\ [^7^(1 — 7s)?] and /I is the renormalization scale. With Wtm 
fermions the renormalization of this operator is multiplicative [16]. The ratio of renormalized 
fi-parameters evaluated in QCD is expected to approach unity as 1 / /I/, — > 0, This follows from 
standard HQET arguments, which also predict the leading deviations for asymptotically small 
l/Zl/j-values to be of the order l/log^/Agco)- Such corrections to the power scaling in 1/p.h 
are expected to be tiny in the /!/, -range of our data and can be estimated in PT by matching HQET 
to QCD. A possible way of doing so is to remove them from the data by considering the modified 
ratio _ 



(4.2) 



where the C-factors ratio contains the info on the l/log( J U/ 1 )-corrections at a fixed order in RG- 
improved PT. We consider here HQET-to-QCD matching only at tree level and LL order in PT 
(thereby avoiding the complications of operator mixing in HQET [17]), and confirm the impact 
of l/log(/I/ ! )-corrections on the final results to be at the level of one standard deviation [18]. In 
the left plot of fig. 4 we show the chiral-continuum limit for the ratio (O s at the heaviest quark 
mass. Both the chiral and continuum extrapolation is well under control. See also ref. [18] for 
additional details. In the right plot we show the heavy-quark mass dependence and the formula 
(O s {Ji h ) = 1 + c\{X) /~pf l + C2W describes the numerical results very well. Applying the ratio 
method one finally gets 

Bf{Ji b ) = 0-90(5), (4.3) 

where the error is the sum in quadrature of statistical and systematic uncertainties. To compute 
BB d , as we have done for the decay constants, we define a ratio for Bb s /Bb cI - Also for this ratio 
the chiral-continuum extrapolation and the heavy-quark mass dependence is very smooth and well 
under control, leading to the preliminar estimate of 

|^ = 1.03(2), < S fe)= 0.87(5). 

Given this result we can also give a preliminary estimate for the parameter 



(4.4) 



fB s jB~B !: 



1.21(6), 



fByBg 

where the error is the sum in quadrature of all uncertainties. 
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Figure 4: Left plot: chiral-continuum extrapolation of the ratio of the bag parameters, co s , for the heaviest 

(7) 

quark mass analyzed, flY '. Right plot: non-perturbative heavy -quark mass dependence of the ratio function 
C0 s (JX h ,X), for X — 1.1784 in the continuum. 
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